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SUMMARY 


Numerous  experimental  results  have  demonstrated  similar  behavior  of  the  velocity 
profile  in  the  near-wall  region  of  the  turbulent  boundary  layer  using  Prandtl's  "Plus" 
scaling  variables.  However,  the  implications  for  similarity  behavior  of  the  near-wall 
turbulent  boundary  velocity  profiles  using  Prandtl's  Plus  scaling  variables  have  not  been 
carefully  explored.  In  the  following  report,  we  apply  the  momentum  balance  type 
approach  to  study  velocity  profile  similarity  using  Prandtl's  Plus  scaling  variables  for  the 
wall-bounded  turbulent  boundary  layer.  It  is  shown  that  the  Plus  scaling  variables  will  in 
fact  satisfy  the  requirements  for  similarity  based  on  this  approach  so  long  as  the  friction 
velocity  values  in  the  flow  direction  are  proportional  to  l/(x -x0)  where  x  is  the  distance 
along  the  wall  in  the  flow  direction  and  x0  is  a  constant.  Experimental  results  are 
examined  and  found  to  confirm  that  certain  datasets  we  tested  do  in  fact  have  the 
friction  velocity  values  behaving  as  l/(x-x0).  For  these  datasets,  we  show  the  Plus 
scaling  variables  satisfy  all  of  the  conditions  for  similarity  using  the  momentum  balance 
type  approach.  However,  the  same  velocity  profile  dataset  plots  also  confirm  that 
Prandtl's  Plus  scaling  variables  do  not  show  whole  profile  similarity  in  general.  Hence, 
we  conclude  that  the  scaling  variables  discovered  by  the  momentum  balance  type 
approach  as  presently  constituted  are  a  necessary  but  not  sufficient  condition  for 
velocity  profile  similarity. 
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1.  INTRODUCTION 


Beginning  with  the  pioneering  work  of  Reynolds  [1],  there  has  been  a  concerted  effort 
to  find  coordinate  scaling  parameters  that  make  the  scaled  velocity  profiles  and  shear- 
stress  profiles  taken  at  different  stations  along  the  wall  in  the  flow  direction  to  appear 
to  be  similar.  For  turbulent  boundary  layers,  the  early  search  for  similarity  was  mostly 
unsuccessful.  This  led  to  the  practice  of  trying  to  find  similarity  in  subregions  of  the 
whole  profile.  Perhaps  the  most  successful  case  of  partial  similarity  was  developed  by 
von  Karman  [2]  and  Prandtl  [3]  and  is  known  as  the  Logarithmic  Law  or  "Log  Law."  The 
Log  Law  states  that  in  a  region  adjacent  to  the  wall  of  the  turbulent  boundary  layer  the 
average  velocity  in  the  streamwise  x-direction  of  a  wall-bounded  turbulent  flow  is  given 
by 


u(x,y) 


=  -In 

K 


i.  fyO 


+  B 


V  v  ) 


(1) 


where  y  is  the  height  perpendicular  to  the  solid  surface  ,  V  is  kinematic  viscosity,  re 
and  B  are  constants,  and  uT  is  the  Prandtl  velocity  scaling  parameter,  the  so-called 


friction  velocity. 


Experimental  plots  of  the  wall-bounded  turbulent  boundary  layer  velocity  profiles  taken 
at  various  stations  along  the  wall  using  uT  as  the  velocity  scaling  parameter  and  v/uT  as 

the  length  scaling  parameter  (the  so-called  "Plus"  scaling  factors)  have  consistently 
demonstrated  similarity-like  behavior  in  the  near-wall  region  of  the  profile.  There  has 
been  much  debate  as  to  whether  this  is  due  to  the  fact  that  the  Log  Law  or,  as  some 
believe,  a  Power  Law  that  is  universally  applicable  for  all  wall-bounded  turbulent  flows. 
In  any  case,  one  area  that  has  not  been  explored  is  the  similarity  implications  using  the 
friction  velocity  as  the  key  scaling  variable.  For  example,  the  question  of  whether  the 
Prandtl's  Plus  scaling  variables  even  satisfy  the  conditions  for  similarity  has  not  been 
studied. 


The  standard  approach  for  studying  similarity  is  to  use  the  dimensionless  momentum 
balance  equation.  The  search  for  similarity  scaling  behavior  for  the  turbulent  boundary 
layer  using  this  approach  began  with  the  experimental  and  theoretical  work  of 
Clauser  [4],  Using  the  friction  velocity  uT  as  the  velocity  scaling  variable  and  the 
displacement  thickness  S1  as  the  length  scaling  variable,  Clauser  predicted  that 


equilibrium  (similar)  boundary  layers  are  only  obtained  for  the  nonzero  pressure 
gradient  case  when  the  quantity 


is  a  constant. 


Pt  =  ~ 


8 1  dp 


(2) 


pu;  dx 

In  this  equation  p  is  the  density  and  pe  is  the  pressure  at  the  boundary 


layer  edge.  Rotta  [5]  and  Townsend  [6]  subsequently  developed  some  additional 
theoretical  considerations  for  turbulent  boundary  layer  similarity  using  this  approach. 
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More  recently,  Castillo  and  George  [7]  used  this  approach  and  found  that  similarity  will 
exist  only  when  the  parameter 

A  =  -  Sdu-  /d*  (3) 

us  dS/dx 

is  a  constant.  In  this  equation  8  is  an  as  yet  unspecified  thickness  scaling  variable  and 
us(x )  is  an  as  yet  unspecified  scaling  velocity.  Following  a  similar  scheme  of  keeping 
the  thickness  scaling  variable  and  the  velocity  scaling  variables  as  unspecified, 
Weyburne  [8]  added  the  dimensionless  Reynolds  stress  transport  equation  to  obtain 
additional  restrictions  with  the  result  that  it  was  shown  that  in  the  outer  region  of  the 
turbulent  boundary  layers  must  be  a  linear  function  of  the  type  8  =  a1(x-x0)  and  us 

must  beus  ccue  and  a  power  law  function  of  the  type  ue  =a2(x-x0)m ,  such  that  al,a2, 
m,  and  x0  are  constants. 

What  is  relevant  here  is  that  in  none  of  these  previous  momentum  balance  type 
approaches  have  explicitly  tested  using  uT  as  the  velocity  scaling  parameter  and  v / uT  as 

the  length  scaling  parameter  to  find  out  whether  they  satisfy  the  conditions  for 
similarity  of  the  velocity  profiles  of  a  turbulent  boundary  layer.  In  what  follows,  we  take 
on  this  task.  Using  the  momentum  balance  equation  type  analysis,  we  first  confirm 
theoretically  that  the  unknown  thickness  scaling  variable  5  and  the  unknown  scaling 
velocity  us(x )  must  behave  in  a  manner  consistent  with  the  sink  flow/wedge  flow 
solution.  It  is  shown  that  the  Plus  scaling  variables  will  in  fact  satisfy  this  requirement 
for  similarity  based  on  the  momentum  balance  equation  approach  so  long  as  the  friction 
velocity  values  in  the  flow  direction  are  proportional  to  l/(x -x0)  where  x  is  the  distance 
along  the  wall  in  the  flow  direction  and  x0  is  a  constant.  Experimental  results  are 
examined  and  found  to  confirm  that  certain  datasets  we  tested  do  in  fact  satisfy  all  the 
conditions  of  whole  profile  similarity  based  on  the  momentum  balance  type  approach 
requirements.  However,  the  same  plots  also  confirm  that  Prandtl's  Plus  scaling 
variables  based  plots  do  not  show  velocity  profile  similarity  over  the  whole  boundary 
layer  region.  Hence,  we  conclude  that  the  scaling  variables  discovered  by  the 
momentum  balance  type  approach  are  a  necessary  but  not  sufficient  condition  for 
velocity  profile  similarity. 


2.  SIMILARITY  EQUATIONS 

The  theoretical  guidance  for  discovering  similarity  scaling  laws  for  the  turbulent 
boundary  layer  started  with  Clauser  [4]  who  looked  at  the  implications  of  similarity 
using  the  x-momentum  balance  equation.  For  a  2-D  incompressible  turbulent  boundary 
layer  that  is  steady  state  on  the  mean,  the  Reynolds-averaged  stream-direction 
component  (x-direction)  of  the  momentum  balance  for  flow  along  a  plate  is  given  by 
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du  du 

, - ^  y - 

ax  ay 


d 

dx  1 


+  —  uu  + 


d 

dy 


1  8P  d2u 

\ uv \  = - +  v — - 

p  dx  dy 


(4) 


where  the  bar  above  a  variable  is  the  Reynolds  average  operator  and  the  tilde  operator 
designates  the  instantaneous  velocity.  What  Clauser  [4]  did  was  to  cast  Eq.  4  into 
reduced  units  in  a  manner  similar  to  what  Blasius  [9]  and  Falkner  and  Skan  [10]  did  for 
laminar  flow  over  a  flat  plate.  What  one  finds  in  doing  this  reduction  is  that  each  term 
in  Eq.  4  can  be  reduced  to  a  product  of  a  x-functional  and  a  ij  -functional  where  ij  is  the 
reduced  y-variable  given  by  rj  =  y/8.  For  similarity  we  must  have  that  the  77 -functional 
terms  do  not  change  as  one  moves  along  the  plate  and  that  the  various  x-functional 
terms  from  Eq.4  change  in  a  proportional  manner  as  one  moves  along  the  plate. 
Equivalently,  we  can  divide  through  by  one  of  the  x-grouping  terms  and  similarity  would 
then  require  that  the  x-grouping  ratios  be  constants.  Clauser  [4]  assumed  that  for  the 
turbulent  boundary  layer  flow  over  a  plate  that  the  length  scaled  as  81  and  the  velocity 

scaled  as  uT  which  then  led  to  Eq.  2  as  one  of  the  x-grouping  ratio  terms. 


A  generalized  process  of  transforming  x-momentum  balance  equation  starts  by  using  5 
and  us  as  the  unknown  similar  scaling  variables  for  the  length  and  velocity.  These 
unknowns  must  be  functions  of  x  but  not  y.  Next  we  assume  that  the  stress  terms  can  be 
separated  into  the  product  of  an  x-dependent  functional  and  an  77 -dependent  functional 
as 

uv  =  uv(x)g12(rj)  and  u2  =  uu(x)  gn(ri)  .  ^ 

The  task  of  casting  x-momentum  balance  and  the  Reynolds  stress  transport  equations 
into  reduced  units  using  these  length,  velocity,  and  stress  scaling  variables  is  detailed  in 
Appendix  1. 


The  results  of  the  transformation  are  that  the  x-momentum  balance  equation  (Eq.  4) 
reduces  to 


A 


— -/'  -T1rjg,11  +  T2gu  +  T3g,n  =  4> 


1- A 


(6) 


where /is  a  dimensionless  function  satisfying  the  Stream  function  (Eq.  1.2),  where  the 
prime  represents  differentiation  with  respect  to  77,  and  where 

Sucd8  .  5  dujdx  ,  u.du.ldx  A 

(p  =  -- 


a  = 


Tl  = 


82  dus 

v  dx 
uu(x) 

<4(1- A) 


8  dujdx 

-  +  — - —  ,  A  = - SJ- — 

v  dx  us  d8/dx 

duu(x)/dx 

us  dujdx  +  i^u]  /S^jd8/dx 


,  and 


r3  = 


us  dujdx  1-A 
uv(x) 

8us  dujdx  +  u2d8/dx 


(7) 


Similarity  of  the  velocity  profiles  means  that  the  parameters  appearing  in  Eq.  6,  such  as 
a,  A,  tx,  t2,  etc.  need  to  be  constant  as  one  moves  along  the  wall  in  the  flow  direction. 


We  note  that  the  x-direction  momentum  balance  equation  is  but  one  of  a  handful  of 
equations  that  govern  the  flow  behavior.  In  an  earlier  report,  Weyburne  [6]  pointed  out 
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that  another  important  flow  governing  equation  appropriate  to  wall-bounded  turbulent 
flow  similarity  was  the  Reynolds  stress  transport  equation  given  by 


duv  duv  ^—du  d2uv  d 

u - h  v - 2  uv - v - - — l - 

dx  dy  dy  dy 2  dy 


r 


uvv  +  - 


v 


Pu 

P 


„  du  dv  P 

+  2v - 

dy  dx  p 


2 du  d\ 7 
v  dy  dx 


=  0 


(8) 


It  is  generally  known  that  turbulent  velocity  profile  similarity  will  only  occur  if  the 

Reynolds  shear  stress  profiles  also  show  similarity.  Hence  it  makes  sense  to  consider 

the  similarity  implications  of  Eq.  8  in  addition  to  the  x-direction  momentum  balance 

equation.  If  one  applies  the  same  reduction  technique  used  above  on  the  momentum 

balance  equation,  one  finds  that  the  Reynolds  shear  stress  equation  becomes 

~Kf  g'u  -£Jlf'gn+T4  f'gu  ~  s  fg'n  +  2gn  f"  +  Zg'n  +  {3  additional  terms }  =  0  ,  (9) 

where  the  three  additional  terms  are  not  written  out  expressly  because  they  cannot  be 

written  in  terms  of  uv  (see  Appendix  1)  and  where 

v  dS  ,  8  du 

Z  =  —  >  e  =  —  ,and/c  = - —  .  (10) 


8uc 


dx 


us  dx 


2.1  Near-Wall  Flow  Similarity 

The  requirements  for  similarity  using  the  momentum  balance-type  approach  consists  of 
ensuring  the  parameters  identified  in  Eqs.  7  and  10  are  constant  as  one  moves  along  the 
wall  in  the  flow  direction.  By  inspection  one  solution  is  to  have 

x  i  \  l  \  (11) 

S  =  a1(x-x0)  ,  us(x)  =  — —  , 

x-x0 

where  alt  a2,  and  x0  are  constants.  This  is  the  so-called  sink  flow  (or  wedge  flow) 
solution. 


For  flows  with  a  pressure  gradient  in  the  flow  direction,  the  momentum  balance 
requires  that  ^  be  a  constant.  Taking  A  as  a  constant,  the  general  solution  to  have  (j) 

(Eq.  7)  equal  to  a  nonzero  constant  is  given  by  us  =  ^a  +  b  u2e  where  o  and  b  are 
constants.  With  only  a  slight  loss  of  generality  (taking  a= 0),  we  see  that  for  similarity, 
this  reduces  to  us  cc  ue  as  first  pointed  out  in  this  context  by  Castillo  and  George  [4], 

2.2  Outer  Region  Flow  Similarity 

In  the  outer  region  of  a  turbulent  boundary  layer  on  a  wall  we  will  assume  that  the 
viscous  forces  are  negligible.  This  means  that  the  viscous  terms  in  Eqs.  4  and  8  are 
negligible  (eliminating  a  and^).  Therefore,  the  functional  form  that  us  and  5  may 
take  is  now  governed  by  the  fact  that  A,  <j) ,  s ,  and  k  from  Eqs.  7  and  10  must  be 
constants.  This  means  that  the  functional  form  that  8  may  take  is  a  linear  function  of 
the  type  S  =  al(x-x0)  such  that  al  and  x0  are  constants.  This  reduces  k  to 


(12) 
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a1(x-x0)  clu 


us  dx 

which  has  the  general  solution  of 

us{x)  =  a2(x-x0 ) 


=  «  , 


where  a,  is  a  constant. 


(13) 


For  flows  with  a  pressure  gradient  in  the  flow  direction,  the  momentum  balance 
requires  that  ^  be  a  constant.  The  general  solution  to  have  ^  (Eq.  7)  equal  to  a  nonzero 

constant  is  given  by  us  =  sja  +  b  u2e  where  a  and  b  are  constants.  With  only  a  slight  loss 
of  generality  (taking  o= 0),  we  see  that  for  similarity,  this  reduces  to  us  ccue  as  first 
pointed  out  in  this  context  by  Castillo  and  George  [4], 

3.  PRANDTL  PLUS  SCALING  SIMILARITY 

The  exact  requirements  for  momentum-balance  type  whole-profile  similarity  can  now 
be  applied  to  the  Prandtl  Plus  friction  velocity-based  scaling  parameters.  Consider  first 
the  near-wall  region.  Substituting  the  length  scale  Sccl/uT  and  the  velocity  scale 
us  ccur  into  Eqs.  7  and  10,  then  it  is  possible  to  show  that  similarity  in  the  near-wall 
region  requires  that  the  friction  velocity  must  be  a  function  of  the  type 

a  (14) 


uT(x)  = 


x-xn 


where  a  and  x0  are  constants.  That  is,  if  the  friction  velocity  behaves  as  Eq.  14,  then  the 
Prandtl  Plus  scaling  parameters  insure  all  of  the  parameters  in  Eqs.  7  and  10  are  equal  to 
a  constant  value.  Next,  consider  the  Outer  region.  In  the  Outer  region,  the  parameters 
given  in  Eqs.  7  and  10  (neglecting  a  and^,  and  setting  k  =  -  1)  will  also  be  constant 
value  so  long  as  the  friction  velocity  follows  Eq.  14.  Hence,  according  to  the  momentum 
balance  type  approach,  all  wall-bounded  turbulent  similarity  requirements  are  satisfied 
using  the  Prandtl  Plus  scaling  variables  so  long  as  the  friction  velocity  behaves  as  Eq.  14 
and  (j)  is  a  constant. 


Having  (j)  constant  will  have  an  additional  constraint  on  the  similarity  requirements  for 
the  Prandtl  Plus  scaling  parameters.  We  note  this  parameter  subsumes  the  Euler 
equation  given  as 

1  dpe  __  due  (15) 


=  ua 


p  dx  dx 

which  therefore  only  applies  to  flows  with  a  pressure  gradient.  For  these  cases,  as 
noted  above,  the  momentum  balance  type  similarity  requires  us  xue .  Since  we  are 
testing  whether  uT  is  also  a  similarity  parameter,  then  we  must  have 
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constant. 


(16) 


the  so-called  Rotta  constraint. 


4.  EXPERIMENTAL  VERIFICATION 

We  are  now  in  a  position  to  answer  the  question  as  to  whether  the  Prandtl  Plus  friction 
velocity-based  scaling  parameters  satisfy  the  momentum-balance  type  whole-profile 
similarity  requirements.  A  number  of  datasets  from  the  literature  were  selected  based 
on  satisfying  Eq.  16.  In  Figs.  la-7a  we  plot  some  sets  of  velocity  profiles  in  which  the 
plot  axes  are  put  in  reduced  units  using  the  Prandtl  Plus  friction  velocity-based  scaling 
parameters.  In  these  figures  the  wall  is  located  at  y+  =  0 .  In  all  cases  there  is  a  region 
near  the  wall  for  which  all  of  the  profiles  overlap,  i.e.  they  display  near-wall  similarity 
behavior.  The  question  we  want  to  answer  is  whether  these  same  datasets  also  satisfy 
the  condition  for  whole-profile  similarity,  that  is  does  the  friction  velocity  values  fit  to 
Eq.  14  and  is  Eq.  16  satisfied.  First,  consider  the  friction  velocity  fits.  In  Figs.  lb-7b  we 
plot  the  experimental  friction  velocity  as  points  (+)  versus  the  station  location  for  each 
corresponding  dataset.  The  red  line  in  each  plot  is  the  least  squares  fit  to  the  points 
using  Eq.  14.  In  general,  the  fits  are  good  indicating  that  Prandtl  Plus  scaling  factors  do 
in  fact  satisfy  the  condition  for  constant  parameter  ratios  (Eqs.  7  and  10). 

For  the  flows  considered  in  Figs.  1-5,  which  do  have  a  pressure  gradient  present  in  the 
flow  direction,  we  must  also  have  Eq.  16  hold  if  uT  is  to  be  a  similarity  parameter.  This 

condition  is  checked  in  Figs.  lc-5c  which  plots  this  velocity  ratio  along  with  average 
velocity  ratio  value  as  the  red  line.  As  is  evident  from  the  plots,  this  constraint  does  in 
fact  hold  for  these  datasets  (recall  that  in  fact  these  datasets  were  selected  based  on 
this  criterion).  It  is  necessary  to  note  that  datasets  presented  in  Figs.  6  and  7  are  for  flat 
plate  flows  for  which  dp/dx=0,  and  therefore  the  (j)  parameter  constraint  is  not 
applicable. 

The  experimental  results  presented  above  indicate  that  all  of  these  datasets  satisfy  the 
conditions  for  similarity  using  the  momentum  balance  type  approach  when  using  ur  as 

the  velocity  scaling  parameter  and  1  / uT  as  the  length  scaling  parameter.  Flowever,  it 

also  apparent  from  Figs.  la-7a  that  these  datasets  do  not  in  general  display  velocity 
profile  similarity  over  the  whole  profile. 
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Figure  la:  The  solid  lines  are  eight  Clauser  [4]  Figure  lb:  Clauser  [4]  friction  velocity  (+)  along 
profiles  plotted  in  Plus  units.  with  the  o/(x-x0)  fitted  line. 
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Figure  lc:  Clauser  [4]  velocity  ratio  (+)  along 
with  the  average  value  as  the  red  line. 
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0  3000 


y 

Figure  2a:  The  solid  lines  are  six  Herring  and 
Norbury  [11]  profiles  plotted  in  Plus  units. 


X  (ft) 


Figure  2b:  Herring  and  Norbury  [11]  friction 
velocity  (+)  and  the  o/(x  -x0)  fitted  line. 


X(ft) 

Figure  2c:  Herring  and  Norbury  [11]  velocity  ratio  (+) 
along  with  the  average  value  as  the  red  line. 
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Krogstad  [12]  profiles  plotted  in  Plus  units.  velocity  (+)  along  with  the  o/(x-x0)  fitted  line. 


x(m) 

Figure  3c:  Skare  and  Krogstad  [12]  velocity  ratio  (+) 
along  with  the  average  value  as  the  red  line. 
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y 

Figure  4a:  The  solid  lines  are  the  Bradshaw 
and  Ferriss  [13]  profiles  plotted  in  Plus  units. 


X  (ft) 

Figure  4b:  Bradshaw  and  Ferriss  [13]  friction 
velocity  (+)  along  with  the  o/(x-x0)  fitted  line 
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Figure  4c:  Bradshaw  and  Ferriss  [13]  velocity  ratio 
(+)  along  with  the  average  value  as  the  red  line. 
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Figure  5a:  The  solid  lines  are  twelve  Jones, 
Marusic,  and  Perry  [14]  profiles  plotted  in  Plus 
units. 


x  (mm) 

Figure  5b:  Jones,  Marusic,  and  Perry  [14]  friction 
velocity  (+)  along  with  the  o/(x-x0)  fitted  line. 


Figure  5c:  Jones,  Marusic,  and  Perry  [14]  velocity 
ratio  (+)  along  with  the  average  value  as  the  red  line. 
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Figure  6a:  The  solid  lines  are  the  Smith  and 
Smits  [14]  profiles  plotted  in  Plus  units. 


Figure  7a:  The  solid  lines  are  the  Wieghardt 
and  Tillmann  [15]  profiles  plotted  in  Plus  units. 


Figure  6b:  Smith  and  Smit  [14]  friction  velocity 
(+)  along  with  the  o/(x-x0)  fitted  line. 


Figure  7b:  Wieghardt  and  Tillmann  [15]  friction 
velocity  (+)  along  with  the  o/(x-x0)  fitted  line. 
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5.  DISCUSSION 


The  implications  of  the  above  results  should  be  viewed  in  the  context  of  two 
observations  which  stand  out:  1)  in  spite  of  satisfying  all  of  the  requirements,  the  set  of 
profiles  given  in  Figs.  la-7a  clearly  demonstrate  that  the  Prandtl  Plus  scaling  parameters 
do  NOT,  in  general,  result  in  similarity  flow  behavior  (except  at  the  near-wall),  and  2) 
while  the  length  scale  Sccl/ur  and  the  velocity  scale  usqcut  have  the  correct  x- 
behavior  (Eq.  14),  it  is  necessary  to  point  out  that  only  two  of  the  geometries  considered 
in  Figs.  1-7  conform  to  the  traditional  sink  flow/wedge  flow  geometry. 

Consider  the  first  observation,  that  is,  the  Prandtl  Plus  scaling  parameters  do  NOT  in 
general  result  in  similarity  flow  behavior  in  spite  of  satisfying  all  of  the  conditions  for 
similarity  using  the  momentum  balance  type  approach.  First  of  all  we  must  point  out 
that  the  datasets  chosen  above  were  not  picked  at  random.  They  were  picked  because 
the  experiments  in  question  were  intentionally  (or  unintentionally)  set  up  specifically  to 
obtain  similar  profiles.  In  fact  only  a  handful  of  wall-bounded  turbulent  boundary  layer 
datasets  fall  into  this  category.  In  any  case,  given  the  examples  above  which  do  satisfy 
all  of  the  constraints  but  do  not  show  whole-profile  similarity,  it  must  be  concluded  that 
the  scaling  variables  discovered  by  the  momentum  balance  type  approach  are  a 
necessary  but  not  sufficient  condition  for  velocity  profile  similarity. 

The  implications  of  this  are  rather  perplexing.  Flow  is  it  possible  to  satisfy  the  conditions 
for  similarity  using  the  momentum  balance  type  approach  and  yet  not  have  similar 
profiles?  One  possibility  is  that  there  is  some  other  factor  that  we  are  not  capturing 
with  the  momentum  balance  type  approach.  Consider  the  results  of  Herring  and 
Norbury  [11]  and  Jones,  Marusic,  and  Perry  [14]  shown  in  Figs.  2  and  5.  It  is  apparent 
from  Fig.  2a  and  Fig.  5a  that  there  are  stations  along  the  flow  direction  that  do  in  fact 
show  similarity-like  behavior  over  the  whole  boundary  layer.  It  is  important  to  note  that 
these  two  datasets  are  the  only  ones  studied  here  that  actually  corresponds  to  the  sink 
flow/wedge  flow  geometry.  Flence  it  is  not  unexpected  that  they  show  similarity 
behavior.  What  is  unexpected  is  that  the  other  datasets  do  not  show  similarity  in  spite 
of  satisfying  all  the  requirements.  Thus  it  may  be  that  there  is  some  other  factor  that 
we  are  not  capturing  with  the  momentum  balance  type  approach  that  explains  why 
these  sink  flow  results  show  similarity  behavior  and  the  other  results  do  not. 

The  second  related  observation  is  also  perplexing  in  that  the  scaling  velocity  under 
consideration  behaves  as  l/(x-x0)  but  not  all  of  the  flow  geometries  are  the 
traditional  sink  flow/wedge  flow  type  geometry.  This  perplexing  result  may  be  due  to  a 
matter  of  definition.  Traditional  sink  flows/wedge  flows  are  associated  with  laminar 
flow  between  two  converging  planes  which  will  have  potential  flows  for  which  the  outer 
velocity  ue  will  behave  as  l/(x-x0).  Since  the  datasets  above  all  obey  the  Rotta 
constraint  and  we  have  shown  that  the  uT  values  in  the  flow  direction  behaves  as 
l/(x-x0),  then  it  follows  that  ue  will  also  behave  as  l/(x-x0),  results  we  have 
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confirmed  but  have  not  shown.  Hence,  not  all  flows  for  which  the  velocity  ue  scales  as 
l/(x-x0)  corresponds  to  the  traditional  sink/wedge-type  flow  geometry. 


6.  CONCLUSION 

Prandtl's  "Plus"  scaling  variables  were  examined  as  velocity  profile  scaling  parameters 
using  the  momentum  balance  type  approach  to  similarity.  It  was  shown  that  certain 
datasets  can  be  found  for  which  the  Plus  scaling  variables  do  in  fact  satisfy  all  of  the 
requirements  for  similarity  based  on  this  approach.  However,  the  same  velocity  profile 
datasets  also  confirm  that  Prandtl's  Plus  scaling  variables  do  not  always  show  whole 
profile  similarity.  Hence,  we  conclude  that  the  scaling  variables  discovered  by  the 
momentum  balance  type  approach  are  a  necessary  but  not  sufficient  condition  for 
velocity  profile  similarity. 
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APPENDIX  1 


In  this  section  the  relevant  equations  for  the  momentum  balance  type  approach  to 
velocity  profile  similarity  are  derived.  Suppose  that  the  x-axis  is  placed  in  the  plane  of 
the  wall  and  that  the  y-axis  is  at  right  angles  to  this  wall  into  the  fluid  layer. 
Furthermore,  suppose  only  steady  state  solutions  are  considered.  We  start  with  the 
Reynolds  decomposition  into  a  mean  (or  average)  component  and  a  fluctuating 
component 

u  =  u  +  u , 

(1-1) 

V  =  v+v, 
p  =  p+p , 

where  the  average  components  are  u,  the  x-velocity,  v,  the  y-velocity,  and  p,  the 
pressure.  The  fluctuating  components  are  u,  v,  and  p  . 

The  equation  for  the  mass  conservation  requires 


du  dv 
dx  dy 


=  0 


du  dv 
dx  dy 


-  0  and 


du  OV 
dx  dy 


=  0. 


For  a  two-dimensional,  incompressible,  turbulent  boundary  layer,  we  ASSUME  the  x- 
component  of  the  momentum  balance  is  given  approximately  by 


du  du  d 

u —  +  v — +  — 
dx  dy  dx 


u2  >  + 


d  <— 


1  dp 

= - -  +  v- 


dx  dy  dx  (  '  dyy  ’  p  dx  dy 

where  p  is  the  density  and  V  is  the  kinematic  viscosity.  Next,  we  introduce  the 
Reynolds  stress  transport  equation  given  by 


duv  duv  „ —  du 

u - h  v - 2  uv 

dx  dy  dy 


_  Pu  du  dv  P  f  du  dv  \  .  M  .. 

uvv  +  —  +2v - +  — [  =  0  .  (1.4) 

p  dy  dx  p\dy  dx 


In  order  to  reduce  these  equations  to  dimensionless  equations  we  start  by  defining  a 
dimensionless  independent  variable 

_  _  y 


where  the  function  S(x)  is  an  as  yet  unknown  boundary  layer  thickness  scaling 
parameter  that  is  only  a  function  of  x.  Furthermore,  we  define  a  stream  function 
y/(x,y)  in  terms  of  a  product  of  functions  as 

i//(x,y)  =  S(x)us(x)f(x,Tj)  ,  (1.6) 

where  f(x,rj )  is  a  dimensionless  function,  and  us(x )  is  the  as  yet  unknown  scaling 
velocity.  We  ASSUME  that  the  stream  function  satisfies  the  conditions 
u  =  dip{xty)  v  =  dip(x,y) 
dy  ’  dx 

This  means  that 
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dx 

where  we  have  used  the  fact  that 

dr/  _  d  [  y 
dx  <3xjr)(x) 

Hence  the  stream-wise  velocity  becomes 


dx  dx 


dy  r, 

u  =  —  =  uj 
dy 

where  the  prime  indicates  differentiation  with  respect  to  77  and  where 

dr/  _  S  j  y  1  _  1 

dy  dy\s(x)\  8 


(1.10) 


(1.11) 


Substituting  the  above  variable  property  similarity  variables  into  the  x-component  of  the 
momentum  equation  (Eq.  1.3),  starting  on  the  left-hand  side,  we  have 


=  ud±fa+ulJ,df 

dx  dx  dx 


The  next  term  is 


Us  diSUs}  ffn  u2  rndf 

8  dx  s  dx 


Combining  these  terms 

[ax  dy  j  sdx  8  dx  5  dx  s  dx 

The  next  step  is  to  transform  the  viscous  component  in  Eq.  1.3  given  by 

=  v%r  . 


(1.12) 


(1.13) 


(1.14) 


(1.15) 


ASSUME  that  the  pressure  in  the  boundary  layer  is  constant  and  equal  to  the  pressure  at 

the  boundary  layer  edge  and  that  the  Euler  equation  given  by 

1  dp  du  (1.16) 

- =  ue  — - 

p  dx  dx 

holds  where  ue{x)  is  the  velocity  at  the  boundary  edge. 


ASSUME  that  the  stress  terms  can  be  separated  into  the  product  of  an  x-dependent 
functional  and  an  77 -dependent  functional  as 

^  =  £12(7)  ,  — —  =  Quid)  ■ 

uv(x)  12V/;  uu(x)  11  W; 

The  Reynolds  stress  terms  reduce  to 


(1.17) 
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(1.18) 


— (T!+— -T 


dx  1  >  dy 


uv\  = 


uu(x)dS„„’  duu(x)„  ^uv(x)„> 


8  dx 


70n 


dx 


0n 


8 


-9 12 


Combining  the  transformed  terms,  we  get  the  x-component  of  the  momentum  balance 
as 


du 


S  r> 2 


dx 


r- 


us  d{Sus 
8  dx 


„df 


fr+uif^-su5r^-+ 


dx 


dx 


(1.19) 


uu 


(x)dS  ,  duu(x)  uv(x )  , 


8  dx 


70n  + 


dx 


„  ,  ,  du  u. 

011+— ^012  =  Ue~T+V^lf 
7  dx  A 


For  similarity,  we  will  ASSUME  that  the  terms  involving  df'/dx  and  df  / dx  are 
negligible.  The  next  step  is  the  key  to  the  problem  of  finding  similarity  restrictions  using 
the  momentum  balance  approach.  Specifically,  it  is  apparent  from  examination  of 
Eq.  1.19  is  that  similar  solutions  will  only  be  possible  if  the  variable  groupings  with  an  x- 
dependence  change  in  the  identical  manner.  Equivalently,  we  can  divide  through  by 
one  of  the  x-groupings  and  similarity  would  then  require  that  the  x-grouping  ratios  be 

u  d{8us\ 

constants.  To  this  end  we  divide  Eq.  1.19  through  by  — — 1 The  transformed 

8  dx 

dimensionless  x-component  of  the  momentum  balance  then  reduces  to 


~  Tl70ii  +  r20ii  +  r3012  =  ^ 


where 


8  dus  8us  dS 

a  = - -  +  — 1 — 

v  dx 


h=  2 


v  dx 
uu(x) 


8  dujdx  ,  updue/dx  A 

A  = - — —  ,</>  =  --  ' 


uc  d8/dx 


h=- 


duu(x)/dx 


,and  r,  = 


ucdujdx  1-A 
uv(x) 

3  8usdujdx  +  u]d8/dx 


(1.20) 


(1.21) 


uJl  +  A)  ‘  usdujdx  +  (u]  js^dd/dx 

The  factors  8 ,  us,  uv{x) ,  and  uu(x )  appearing  in  Eq.  1.20  are  as  yet  unknown 
parameters. 


In  a  similar  fashion  we  can  transform  the  Reynolds  Stress  Transport  equation  given  by 

Eq.  1.4.  Substituting  the  above  similarity  variables  into  the  Reynolds  stress  transport 

equation,  starting  on  the  left-hand  side,  we  have 

duv  d{uv(x)g12  (77)}  ,\usdd  <  duv(x)  U-22) 

u—  =  “sf  - - - ~  =  -uv(x)^—?ifg12+u— — fg12  . 

ox  dx  8  dx  dx 


The  next  term  is 

duv  ,  1  d{8us}  df  , 

v =  -w(x)-  1  ’  fg12-uv(x)u—g12 

dy  8  dx  dx 

The  first  term  on  the  right-hand  side,  the  Production  term,  is 
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(1-23) 


(1.24) 


(1.25) 


du  usu  i/(x) 

2  uv —  =  2  — — — g„/ 

dy  8  12 


The  next  step  is  to  transform  the  viscous-like  component  in  Eq.  1.4  given  by 

r?2uv  uv(x) 

dy 2  £2  12 

Combining  the  transformed  terms,  we  get  the  stress  balance  as 

/  ,du  ,  ,ud8  ,  duv(x)  u  dS  ,,  , 

-uv(x)—±fg12-uv(x)^—fg12  +  us  w/  g12 -uv(x)-^—>ifg12  + 
dx  8  dx  ox  8  dx 


(1.26) 


ucuv(x)  „  uv(x)  „  ,  , 

2 - 9u  +  v — hr~Qi2  +  { 3  additional  terms)  =  0  , 

8  8 

where  the  three  additional  terms  are  not  written  out  expressly  because  they  cannot  be 
written  in  terms  of  uv  .  This  is  not  to  say  that  these  terms  can  be  neglected.  In  fact  the 
opposite  is  true;  the  three  additional  terms  in  Eq.  1.26  include  the  numerically 
significant  energy  dissipation  rate  and  the  velocity-pressure  gradient  terms.  However, 
for  the  purposes  of  obtaining  similarity  scaling  information  for  8  and  us,  these  three 
additional  terms  would  require  additional  assumptions  and  the  identification  of 
additional  unknown  parameters.  Hence,  for  the  purposes  herein,  we  will  leave  the 
terms  unused. 


The  last  step  is  to  insure  that  any  of  the  variable  groupings  with  an  x-dependence  are 
constant.  In  order  to  make  the  stress  equation  dimensionless  we  divide  through  by 
umv(x) 

-  to  get  the  transformed  Reynolds  stress  balance  as 

-Kfg[2  -£r/f'g'12+T4  f'Qu  ~s  fg  12  +26,i2  f"  +  X9i2  +  {3  additional  terms}  =  0  .  (1-27) 

where 


v  d8  .  8  du 

y  =  -  ,  e  =  —  ,  and  k  = - 5- 

Su  dx  u  dx 


(1.28) 
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